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1. Introduction 

Large interest has been recently raised by the analytic solution of the 7— unstable 
collective Bohr hamiltonian in the case of a square well potential in the (3 variable 
PP. This situation describes the shape phase transition between spherical and 
7— unstable nuclei and it has been associated with the E(5) group. Similarly an 
approximate separation of variables has been proposed [2] as a solution for the 
transition between spherical and axially deformed shapes which has been referred 
to as X(5). Experimental confirmation of the actual occurrence of such situations was 
found for instance in the examination of the level scheme of 134 Ba for E(5) and 152 Sm 
for X(5) [HIE]. 

The occurrence of dynamic symmetries is associated with systems where the 
hamiltonian is solely written in terms of Casimir operators of a chain of subalgebras. 
The wider algebra of this chain is the one that settles the nature of the problem 
and its general form bears all the information on the system in a compact way. The 
non-unique way in which the algebra is split in chains of subalgebras exhibits some 
peculiar feature of the system leading to exact analytical results, that may help in 
elucidating experimental findings and data trends. Concerning ourselves to nuclei, an 
example of occurrence of dynamical symmetries is provided by the Interacting Boson 
Model. With particular choices of the general SU(6) hamiltonian one obtains limiting 
situations that are analytically solvable (0(6),U(5) and SU(3) cases). Moving around 
in the nuclear chart or along a chain of isotopes, nuclei can undergo phase transitions 
from one dynamical symmetry to another. The interest of E(5) and X(5) lies in the fact 
that they are related to new analytically solvable situations corresponding precisely 
to the critical phase transition points. 

In the specific case of the transition from spherical to 7— unstable nuclei, the E(5) 
description assumes, for the transition potential in the (3 variable, an infinite square- 
well potential [Q, a case that has been generalized by Caprio with the introduction of 
a finite square well 0. Other choices are possible, like the Davidson potential studied 
by Elliott and later by Rowe jZj , that also generates an analytic vibration-rotation 
spectrum. The work of Rowe is very much akin to the one discussed here, displaying 
the same SO(2,l)xSO(5) algebraic structure. 

The purpose of this investigation is to show that, within the condition of 
7— instability, there are other classes of analytically solvable potentials. 

The Bohr hamiltonian [Hj will be our starting point: 



H ~2B 



P 4 dp 1 dp p 2 sin 3 7 d 1 m 7 d 1 
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V{P,l). (1) 



Following the standard procedure, when the potential depends only upon /3, i.e. 
V(P, 7) = U((3), one can separate variables as in Wilets and Jean obtaining a 
system of two differential equations, one containing the 7 variable and the three Eulcr 
angles, the other containing only the /? variable. The spectrum is determined by the 
solution of the latter, namely 
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U(J3)}f(j3)=Ef(J3) (2) 
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where A 2 is the Casimir operator of SO(5) (r being the associated quantum number) 
|Sj and we rewrite here the differential equation in the (3 variable in the second order 
standard form as: 

X "(/3) + {e- u(fi) (T + p l /2)2 + ±}m = (3) 

where x(P) = /? 2 /(/3), while e = t^t-E and u = W-U are the reduced energies 
and potentials. In this note we will point out that this equation displays analytic 
solutions with the choice of the Coulomb and the Kratzer potentials in the quadrupolc 
deformation parameter (3. The former has a very simple form diverging in zero, while 
the latter has been widely used in the early stages of quantum theory to describe 
interactions within ions in configuration space and has a minimum for a finite value 
of p. In both cases we have in fact the possibility to regain, with some simple 
mathematical steps, the well-known Whittaker's standard form for eq. (J2J, and hence 
to obtain analytic solutions. 



2. Coulomb-like potential 

Inserting the potential 

ucifi) = -j, A>0 (4) 

in equation {3} and with the substitutions e = — e, x = 2y/ej3, k = Aj{2%fe) and 
fj, = t + 3/2 the previous equation takes the Whittaker's standard form 

The solution for negative energies, that is regular in the origin, may be found (as in 
HOI El) to be the Whittaker's function M k ^(x) : 

XkA x ) = K,^ +1 h~ x / 2 1 F 1 (n+ l/2-k,2fi + l;x) (6) 
The normalization Af T £ of these states can be obtained by quadrature § from: 

xWfdp = / P 4 f(f3) 2 dp = 1 (7) 
la Jo 

where the volume element in the (3 variable is f3 A df3 [S]. The hypergeometric function 
for x — > 00 is in general proportional to e x and hence diverges. However, it is not 
divergent when it becomes an associated Laguerre polynomial, that is to say when the 
first parameter /.t+1/2 — fc = r + 2 — Aj i2^fe) is a negative integer, — £. This leads 
to a condition that fixes the spectrum as 

A 2 / 4 
(rTi+2) 2 

Now r + £ works as a single quantum number n for the energies, but the shape of the 
wavefunctions depend on r and £ separately. This spectrum is depicted in figure ^ 
where the energy of the first two states have been fixed respectively to and 1 and 
this is sufficient to settle the energy scale e'. The (4 + , 2 + ) doublet with n = 2 has an 

§ The general integral of two hypergeometric from to 00 is not known, but one can find 
the normalization constants in an analytical way since the hypergeometrics reduce to Laguerre 
polynomials. 



^ = 0,2 (8) 
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Figure 1. Spectrum of the Coulomb-like potential. The energy scale (e') is 
chosen by fixing the energy of the first two states respectively to and 1. The 
transition rate for 2q" x — > 0q~ has been fixed to 100. Some selected quadrupolc 
transitions are shown in the figure for simplicity. 



energy of 1.35 that is smaller than the corresponding values for other well-known cases. 
Furthermore this spectrum displays another interesting feature: there is a threshold 
at 1.8 that corresponds to an infinite quantum number. Again this value is smaller 
that the energy of the two-phonon state of the vibrator. 

For the sake of completeness we give the wavefunctions with the reduction of the 
hypergeometric to associated Laguerre polynomials: 

x(x) t£ =K* x2 " +1 ^ x/ \ 2fi + lk L( t ) ( x ) ( 9 ) 
where the denominator is a Pochhammcr symbol. 

3. Kratzer-like potential 

We move now to the study of the Kratzer potential (sec figure |5J) 

„(«__„,(*-»«) (JO, 

where T> represents the depth of the minimum, located in (3q. Wc may also write it as 
uk{P) = uc{P) + B/P 2 for later purpose. Inserting the cited potential in ®, with 
the substitutions e = -e, x = 2^/ef3, k = Vf3 /(^/e) and /i 2 = (r + 3/2) 2 + X>/3 2 , the 
equation takes again the Whittaker's standard form. The regular solutions are again 
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Figure 2. Kratzcr-likc potential. For purpose of illustration the horizontal axis 
has been set in the adimensional variable /3o//3, while the vertical scale is in units 
of the depth of the pocket. The minimum corresponds thus to /3q. 



Whittaker's functions with the new substitutions and the same arguments apply for 
the properties of convergence. Now fi + l/2-k = v /(t + 3/2) 2 + V@l + 1/2- V/3 a /^e 
must be a negative integer, — £. The spectrum worked out from the last requirement 
reads: 

v 2 a 2 



(A + 2 " (\/(t + 3/2) 2 + Vf3 2 + 1/2 + 5 

A 2 1 4 



(11) 



(v/(r + 3/2) 2 + B + 1/2 + 2 

with £ = 0,1,2,.. . The proper set of quantum numbers that characterizes the 
eigenvalues is r, £, L, M. Notice that each t, £ state may be degenerate with respect to 
the angular momentum, according to the Wilcts and Jean rules jHj. The r quantum 
number is contained in A = yir-\- 3/2) 2 + P/3 2 + 1/2 and the £ quantum number has 
the same meaning as in the paper of Iachcllo (shifted by one unity) being connected 
with the zeros of the Whittaker's function that are determined by the zeros of the 
hypergeometric functions and give rise to the different bands. If one displays the 
spectrum given in formula 111 l[l imposing that the ground state (r = 0, £ = 0) is at 
zero energy and that the (r = 1,£ = 0) state has e = 1, it is evident that one has to 
play with the position of the minimum (3q and the depth T> of the potential. In fig. 
El we study the dependence upon the depth of the potential well. It is seen from the 
lowest few states of the first two bands that there are clearly two limiting cases: when 
the depth goes to zero the spectrum becomes equivalent to the already discussed 1/(3 
case, that we wish to call the Coulomb-like limit, while when T> tends to infinity all 
the £ > bands escape to an infinite energy and the £ = band has a spectrum that 
follows a simple /3-rigid, 7-soft rule: e T = t(t + 3)/4. In fig. 0]we study instead the 
evolution of the spectra with /3o for a fixed value of D. Again the same two limits 
are seen when the deformation parameter is very small or very large. There is a 
qualitative equivalence between the case with small /3o, the case with small T> and the 
Coulomb-like case. Moreover when (3q or T> tend to infinity the situation is equivalent. 
Apart from the limiting cases the most interesting situations are realized for spectra 
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Figure 3. Evolution of spectra with fixed /3rj = 0.5 and increasing 27. The 
first two bands (£ = 0,1) are displayed with their lowest states (r = 0,1,2,..). 
The various substates are not displayed for the sake of simplicity. The transition 
showed are 4^" — » 2^~ and 2^ — > 0q lt normalized with B(E2; 2^ — » 0j ) = 
100. 



that have a reasonable (3o (that is usually smaller than 0.5) and a 27 freely adjustable 
that can widen the relative distances in the £ = band as well as move the lowest 
state of the second £ = 1 band. 



4. Transitions 

Electromagnetic transition rates (see ^OHE]) are defined as reduced matrix elements 
of the transition operator T(E\) 

| (^ n ,L t \\T(E\)M f ,T f ,L f ) * 



B{EX;i i ,T i ,L i -> £f,Tf,L f ) (2L +1) 

and the quadrupole transition operator, to the first order, reads: 



T(E2, fx) oc /3 



n( 2 ) 



(2) 
2 



(2) xSinj 



(12) 



(13) 



The matrix elements have been calculated for a few selected transitions and are 
displayed directly on the figs. 11131 and 0] The (3 part of the integration factorizes 
out and the {7, Oi] part has been calculated in the standard way. Taking into account 
all the transitions between a given r state and one of his neighbouring (for instance 
t — 1) states, it suffices to calculate explicitly one transition, the others being obtained 
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Figure 4. Evolution of spectra with fixed 15 = 10 and increasing @q. The first two 
bands (£ = 0, 1) are displayed with their lowest states (r = 0, 1, 2, ..). The various 
substates are not displayed for the sake of simplicity. The selected transition rates 
are the same than in the preceding figure. 



by fixed geometrical factors coming from the 7-angular part |lri| . It is interesting to 
note that the transition from the first 2+ state of the second band to the ground state 
of the system is very small if compared to the reference transition, evidencing almost 
no overlap between the two wavefunctions in the (3 variable. 

For simplicity in the last two figures we have given only the following transition 
rates: 6^ — > 4^~ , 4^ — > 2^~ and 2^ — > 0q 1; relatively to the transition 2^~ — > 0(| 
that has been given the reference value of 100. It is interesting to note that in both 
cases (fixing one parameter and changing the other) the two evidenced limits have also 
common values for B(E2)'s and the two Coulomb-like limits for the Kratzer potential 
show the same values of the exact Coulombdike case in fig. ^ 



5. Algebraic structure 



Wc discuss now the group theoretical interpretation of the hamiltonians that we have 
solved directly in the previous sections. We will consider the quantization of the 
collective model, restricted to quadrupole deformations only (A = 2), by defining pairs 
of canonically conjugate operators on the Hilbert space = and = — i-^- 
(having dropped any h and using a covariant / contravariant notation) . These operators 
obey the Heiscnberg-Weyl commutation relations [a^, n v ] = i5^^ ^1]. If a is defined 
as the vector whose five components arc the operators defined above, the scalar product 
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where the last equivalence is a consequence of the property: a M = a* = ( — 

Notice that 1 = a • a. We need also to consider the parameters et M and p^ which 
play the role of the a M and 7r M in the intrinsic frame of reference, and have the useful 
property that a M = a_ M . Since the transformations between these two sets of variables 
are unitary, the commutation relations are preserved as well as the scalar products 
(for example tt"^ = E^P/J- 

Taking the reduced quantities, equation may be recast in the form 

7T 2 +u(p) -el /(/?) =0 (15) 



n 



where (see [5]) 



1 d ni d A 



2 



- 2 - = p-p= -J^% + W (16) 

For both the Coulomb-like and Kratzer-like potentials the hamiltonian Tt is 
SO(5) invariant. In addition it displays a spectrum generating algebra since one may 
construct three operators [7] that are infinitesimal generators for the corresponding 
group: 

Zi = 4/?(p 2 + ^) Z 2 = p Z 3 = 2(S-f-i) (17) 

which have the commutation relations of the four non-compact isomorphic Lie algebras 
su(l,l) ~ so(2,l) ~ sl(2,R)~ sp(2,R) (using Wybourne's notation for symplectic group 
dimensions [T5]1: 

[Z 1 ,Z 2 ] = -4iZ 3 [Z 3 ,Z 2 ] = -2iZ 2 [Z 3 ,Z 1 ]=2iZ 1 . (18) 

With the potential u([3) = uk{P) (that contains also the Coulomb-like case, when 
B = 0) the operator (3Ti is in fact expressible as a linear combination of the elements 
of the algebra of su(l, 1), namely in the form 

pH = Zx/A-A. (19) 

Now one can define new operators Xi with i — 1,2,3 by means of a linear 
transformation: 

Xi = \{Z X -Z 2 ) X 2 = X -Z 3 X 3 = i{Zi+Zi) (20) 
which satisfy the following commutation relations: 

[X 1 ,X 2 ] = -iX 3 [X 2 ,X 3 ]=iX 1 [X 3 ,X 1 ]=iX 2 . (21) 

The eigenvalue equation for the Bohr hamiltonian is now, having multiplied by (3 on 
the left, 

P(H-e)V = 0. (22) 

Using equation (|19|) and the definitions of the Xi and Zi operators, we rewrite the 
last equation as 

\l + Ae)X 1 -2/4 + (l-4e)X 3 j* = (23) 
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and following the procedure in we can perform a (1,3) hyperbolic rotation of an 
angle 8 to diagonalize the eigenvalue equation and choosing tgh{0) = — (l+4e)/(l— 4e) 
(valid for e < 0) we obtain 

A 

X 3 ^ = ^=*, (24) 



where ^ is the rotated wavefunction. The Casimir operator of the so(2,l) algebra is 
evaluated to be: 

C 2 = A 2 + X- - X+ + B + 2 (25) 

with eigenvalue t(t + 3) + B + 2. The two last equations must be compared with the 
two following eigenvalue equations (for unitary representations D + |17|): 

x 3 1 0,0 = (f-0) | 0,0 

Ca | 0,0 = 0(0 + 1) |0,O- (26) 
This comparison yields a spectrum of the form: 

A 2 /4 

6r,£ = 7 = (27) 

' 4 (V(r + 3/2)2 + 5 + 1/2 + C) 2 

that coincides with the one found from the direct solution of the differential equation 
with a Kratzer-like potential. It also contains as a special case (B = 0) the Coulomb- 
like case. 

The algebra associated with the SO(5) group is the so-called degeneracy algebra 
(according to the definitions in JH|) and the group SO(2,l) is associated with the 
spectrum generating algebra. The relevant chain of subalgebras that gives the labels 
of the set of orthonormal states {| ^raLM}} is explicitly given as |71 114): 

£17(1,1) x SO{5) D U(l) x SO{3) D 50(2) 

A r a f L M [ 6) 

where A is an SU(1,1) lowest weight and a indexes the SO(3) multiplicity. These basis 
diagonalize the hamiltonian 119() given above. We can thus state that the problem 
studied so far displays a SO(2,l)xSO(5) dynamical algebra. 

6. Final remarks 

In this paper we have solved the Bohr hamiltonian for two specific 7— unstable 
potentials that yield analytical solutions. We have given the corresponding spectra and 
wavefunctions in closed form and the most important transition probabilities. Rather 
interesting looks the case of the Kratzer-like potential that may be given in terms of 
two unrelated parameters A and B, as — A/ (3 + B / 1 . Changing the value of B from 
zero to a finite value, one can describe situations ranging from spherical to quadrupolc 
deformed shape. The critical point is actually the potential with B — > 0, that is the 
Coulomb-like limit discussed above. The same dynamical symmetry discussed briefly 
in the former section, is thus effective for a class of different potentials, that, depending 
on some parameter, may describe very different situations. We would like to remark 
that these potentials parallel the case of the family of potentials of the form: 

u H (P)=Ap 2 (29) 

and 

uu08)=«ff(/3) + 5 (30) 



Analytically solvable potentials for "/-unstable nuclei. 



10 



where H stays for harmonic and D for Davidson, that also lead to solvable Bohr's 
hamiltonians and are furthermore both characterized by a SO(2,l)xSO(5) dynamical 
group 0. Other interesting possibilities, such as linear combinations of powers (to be 
solved with Frobenius method), arise in the same spirit of this paper. 

The (3— part of the spectrum of the potentials that have been discussed here was 
combined with the condition of 7— instability. We will show in a forthcoming paper 
that approximate solutions can also be obtained for the same functional dependence on 
(3, but with a 7— dependence favouring axial symmetry. This case will be therefore the 
homologous of the situation associated with the occurrence of X(5) in the approximate 
solution of reference [21 - 
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